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Abstract 



We calculate the axial form factor in the chiral quark soliton (semibosonized 
Nambu - Jona-Lasinio) model using the semiclassical quantization scheme 
in the next to leading order in angular velocity. The obtained axial form 
factor is in a good absolute (without additional scaling) agreement with the 
experimental data. Both the value at the origin and the g-dependence of the 
form factor as well as the axial m.s. radius are fairly well reproduced. 
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Recently, including the 1/N C rotational corrections (next to leading order) in the semi- 
classical quantization scheme, a natural solution for the problem of strong underestimation 
of the axial-vector coupling constant in the leading order in the semibosonized Nambu - 
Jona-Lasinio (chiral quark soliton) model has been found . These corrections lead 

to an enhancement || of order (iV c + 2)/N c for g^ and improve considerably the agreement 
with the experiment. Apparently one should expect that the 1/N C corrections will not be 
small also for the axial form factor and it can change significantly the existing leading or- 
der results j5|. It should be noted that the latter suffer from the fact that g A = Ga(0) is 
strongly underestimated which does not allow for an absolute agreement with experiment. 
It is usually assumed || that the axial form factor is scaled in the same way like the axial 
vector coupling constant gA- This is, however, not necessarily true since the axial form 
factor contains a second contribution, which does not contribute to g A = G A (0). Also in the 
calculation || a low constituent mass M = 363 MeV is used, which differs from the value 
M = 420 MeV extracted from a overall fit [ITT] to the static properties as well as to the 



electromagnetic form factors of the nucleon in the model. Therefore, it is worth to evaluate 
the axial form factor in the chiral quark soliton model using the semiclassical quantization 
scheme with the next to leading 1/N C rotational corrections included, which is the aim of 
the present work. 

We start with the general decomposition of the matrix element of the axial current 
A% = \t r '7 Q 7 M 7 5 : ^ L1 J/ in terms of the corresponding form factors: 

(N(p', O |* f 7 V7fe£*l# (P, 0> = O \GA(q 2 )r + ^rP-q* 1 } 7te£«(p, , (i) 
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where q = p' —p and £ stands for spin and isospin. After some standard manipulations 
we can express the axial form factor as 



G A (q 2 ) = 3 ^ Jd 3 x e «*» [{p t \Al(x)\p T) - ^(p T \A\(x)\p T>], (2) 



where and E are the nucleon mass and energy E = + respectively, and 

\p |) is a proton state of spin up. It should be noted that eq.(|2|) contains matrix elements of 
the space components of the axial current A\, and apparently only the first term in eq.(^) 
contributes to g A - 

For the evaluation of the matrix element of the axial current we use the simplest Sub- 
version of the chiral quark soliton model with up and down quarks, degenerated in mass. It 
is based on a semibosonized Nambu Jona-Lasinio lagrangean J7],[8| : 

C = ^(-tYd tl + m + MU^)ty } (3) 

which includes auxiliary meson fields 

U{x) = e if -*&/U j (4) 

constrained on the chiral circle. The model is non-renormalizable and one needs a cutoff to 
make it finite. This cutoff A and the current quark mass mo are treated as parameters of the 
model and are both fixed in the mesonic sector to reproduce the physical pion mass m n and 
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the pion decay constant /„-. The third model parameter, the constituent quark mass M, can 
be related to the empirical value of the quark condensate but it still leaves a broad range 
for M. Actually, in order to obtain an overall good description of the baryonic properties a 



value around 420 MeV has to be used 11 



In the model the baryons appear as a bound state of N c (number of colors) valence 
quarks coupled to the polarized Dirac sea. Since the model lacks confinement the proper 
way to describe the nucleon is to consider || a correlation function of two iV c -quark currents 
with nucleon quantum numbers at large euclidean time-separation. Accordingly, the nucleon 
matrix element of the axial current \l/^7 7 /i 75 is represented by an euclidean functional 
integral with lagrangean (|3|): 

(N'J\^(0)~f°^ 5 -V{0)\N,p) = lim - / d 3 xd 3 ye~^ + ^ 

2 T^oo Z J 

x JvU Jvm J v¥J n ,(x',T/2) ^ t (0)7 7' 1 75^(0)4(f,-T/2)e-/ d42 * tD ( c/ )* . (5) 



Current is a composite N c quark operator |I(J with nucleon quantum numbers JJ3, TT 3 : 



MZt) = ^ e^-^r^-^ J Js,TT 3 y^frt) . ..* fiHcflfc (2,t) . (6) 

For the evaluation of the path integral in eq.fl^) we follow the line of ref. [0,|lTJ. Here 
we will only sketch the derivation. Integrating out the quarks in it is easy to see that 
the result is naturally split in valence and sea parts. After that we integrate over the meson 
fields U in saddle point approximation - large N c limit. It leads to a stationary localized 
meson configuration (soliton) of hedgehog structure 

U(x) = e i? -* p ( x ), (7) 

which minimizes the effective action 

Tr log D(U) = Tr log[«9 r + h(U)] , (8) 

where the one-particle hamiltonian h is given by 

a ■ V 

h(U) = —— + PMXT* + pm . (9) 

Since the hedgehog soliton field configuration U(x) does not preserve the spin and isospin, as 
a next step we make use of the rotational zero modes to quantize it @. It is done assuming 
a rotating meson hedgehog fields of the form 

U(x, t) = R(t) U(x) R\t) (10) 

with R(t) being a time-dependent rotation SU(2) matrix in the isospin space. It means that 
even in leading order in iV c one has to go beyond the saddle point approximation extending 
the path integral in eq.([|) over all fields of the form (0) - a path integral over R. It is easy 
to see that for such an ansatz one can transform the effective action 

Trlog D(U) = Trlog(D(U)+iQ) (11) 
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in order to separate the angular velocity matrix: 

n = -iR\r)R{r) = l -n a T a . (12) 

The dot stands for the derivative with respect to the euclidean time r. Similarly, the quark 
propagator in the background meson field U can be rewritten as 

= (x^j^^R^W). (13) 

Since the angular velocity is quantized according to the canonical quantization rule, it 
appears as Q a ~ This allows one to consider Q as perturbation and to evaluate any 
observable as a perturbation series in fl which is actually an expansion in 4-. Actually, we 
make essentially use of the expansions 

N c Tr\og[D(U) + ifi] = N c Tr \og[D(U)] + Q sea J drQ 2 a + ... , (14) 

and 

1 = J: + (15) 

D{U)+iQ D{U) D(U) D(U) 1 ' 

up to terms quadratic in Q to arrive at a functional integral over the time dependent ori- 
entation matrices R{r) with an action quadratic in angular velocities |§. In large N c limit, 
the latter corresponds to the hamiltonian of the quantum spherical rotator: 

H rot = J a 2 /(20) , (16) 

where J a is the spin operator of the nucleon and 

6 = Q sea + Q val ~ N c (17) 



is the total moment of inertia |§[12]> including both the valence and the sea quark contribu- 



tions. It means that despite of the fact that in general the path integral over R runs over all 
possible trajectories, in large N c limit the main contribution comes from trajectories close to 
those of the quantum rotator with hamiltonian (0). According to eq. (16) the quantization 
rule (in euclidean space-time) is given by 

(18) 

Due to the collective path integral over R the order of the not-commuting collective operators 
n a (R(r)) = -iTr(R)(r)R(r)r a ) and D ab (R(r)) = ^TT(R'(r)r a R(r)r b ) , (19) 



is strictly fixed by the time ordering^. For given spin J, J 3 and isospin T, T 3 the spin-flavor 
structure of the nucleonic solution can be expressed through the Wigner D function 



1 The details of this procedure can be found in ref. [11] 
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\N,T 3 J 3 )(R) 



■1) T+T *VW+1D T _=J(R) 



(20) 



In the above scheme, the matrix element of the space components of the axial current 
A% includes leading order terms ~ Q° as well as next to leading order ones ~ Q (1/ iV c ) . In 
Minkowski space-time it has the following structure: 



(N\A k 3 (x)\N) = -id ($L(f)7V75r 3 ^(^)) - E^A°( e n) K^/tS^ 

I n 

%i al (x)j°j k j 5 T b $ n (x)) (n \t c \ val) 



l x 



E sign(e n ) 



(21) 



- ^ ftf (e n , e m ) ($t rt (x) 7 V75^$„(f)) (n |r c | m) 



Here $ n and e„ are the eigenfunctions and the eigenvalues of the hamiltonian (|9]). In eq . (pT|) 
first and third terms are valence quark contributions in leading and next to leading order in 
angular velocity, respectively. They are finite and do not need any regularization. The other 
two terms represent the divergent Dirac sea part and need regularization. The regularization 
functions Tlf, Tlf can be found in refs. 



Inserting the result fl2~l|) in eq.(|2|), after some straightforward calculations we obtain for 
the axial factor 



G A (q 2 



M 



N 



E 



r 2 dr \j (qr)A (r) - 32{qr)A 2 {r)\ 



(22) 



where both densities Aq(t) and A2(r), split in valence and sea parts, contain leading and 
next to leading order terms in angular velocity Q: 



3V3 



$ljr)\\T^\\$, m ,(n 



-0(2) \\^val\ 



+i E pK n + 1 <(e n ) ($t (r) WT^ ||$ B (n)) }, (23) 

n=all ) 



and 



Y^l E (^LWI|T %||$ n (r))M||r«||n> 
1 9V2 [ n ^ val e„ - e vai 

-\ E ^(^OKCOIIlj&ll^mCr^nllrWlIm)}. (24) 



n,m=all 



The tensors T^ 1 ' are defined as 



IT 



C « 



(1) ^^.(1)1(0(1)) 



(25) 



and 



(26) 
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The solitonic solution U is found by solving numerically the corresponding equations of 
motion in an iterative self-consistent procedure ||. To this end we use the method of Ripka 
and Kahana [I3J for solving the eigenvalue problem in a finite quasi-discrete basis. 



We also calculate the axial m.s. radius (r ) A given by 



(r 2 ) A =r 6 dG ^ 2 



GU(0) dg 2 



2 



r 4 dr [A (r) + -A 2 (r)] + — T , (27) 



?2=0 G A (0) J ^ 4M^ 

and for convenience, parameterize also the calculated form factors using a dipole fit: 

G A (q 2 ) 
G A (0) 




(28) 



Our results for the axial form factor are displayed on fig.[3] for four different values of 
the constituent mass M, namely 360, 400, 420 and 440 MeV. Similar to g A [0 the main 
contribution to the axial form factor in the leading as well as in the next to leading order 
in angular velocity comes from the valence quarks and the Dirac sea contribution is almost 
negligible. Table [I] contains the corresponding values for the dipole mass and m.s. radius as 
well as for the axial vector coupling constant. Similar to the other nucleon properties |11] a 
value for constituent mass around 420 MeV is preferred. For this mass value the theoretical 
curve agrees fairly well with the experimental data without any scaling. The lowest value of 
360 MeV is almost excluded. On fig.[l] we also show the leading order results. As can be seen 
the two groups of curves differ in both magnitude and slope. Indeed, the extracted dipole 
masses for the axial form factors, calculated in leading order, are larger. In particular, for 
M = 360 MeV the dipole mass in the leading order still agrees with the experiment, whereas 
in the next to leading order it is out of the experimental window. It means that a simple 
scaling behavior of the axial form factor as it is assumed in ref. is not valid. On fig.|2| we 
present the separated A Q - and /^-contributions to the axial form factor. As can be seen the 
A part completely dominates the axial form factor at small momentum transfer. 

In Table [I] we also present our results for the axial m.s. radius calculated from eq.(f2~7|) 
as well as from the dipole fit (in brackets) in comparison with the estimate from the ex- 
perimental dipole fit. As can be seen the values directly calculated from the corresponding 
densities A Q and A 2 are very close to those of the dipole fit which is an indication that even 
at small q 2 the dipole fit is a good approximation to the theoretical curves. 

To summarize, we evaluate the axial form factor in the chiral quark soliton model tak- 
ing into account the next to leading order (1/N C ) rotational corrections in the semi classical 
quantization scheme. These corrections provide not only the enhancement, needed to repro- 
duce the experimental value of the axial vector coupling constant, but also change the slope 
of the axial form factor allowing for an almost perfect absolute (without additional scaling) 
description of the experimental data. The axial m.s. radii are also well reproduced. 
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FIGURES 



FIG. 1. Axial form factor in leading and next to leading order in angular velocity in comparison 
with the experimental data [p. IS]. 



FIG. 2. Separated Aq- and ^-contributions to the axial form factor. 
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TABLES 



TABLE I. Axial properties of the nucleoli, calculated in the NJL model for four different values 
of the constituent mass M = 360, 400, 420 and 440 MeV, compared with experimental values. The 
dipole masses for the form factors calculated in both leading and next to leading order are given. The 
m.s. radius obtained from the dipole fit (in brackets) is also presented. 



Constituent quark mass M [MeV] 

360 400 420 440 exp 
[GeV] 0.910 ± 0.001 1.004 ± 0.001 1.040 ±0.001 1.072 ±0.001 

M%° +Ql [GeV] 0.849 ±0.001 0.956 ± 0.001 0.995 ± 0.001 1.031 ± 0.001 1.05±g;i§ 

< r 2 > A [fm 2 ] 0.70 (0.65) 0.51 (0.51) 0.46 (0.47) 0.43 (0.44) 0.42±g;J| 

g A 1.31 1.24 1.21 1.18 1.26 
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